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ABSTRACT 

Capture cross sections are calculated for ion-polar molecule colli- 
sions of interest in mass spectrometry. These cross section predictions 
are in relatively good agreement with experimental data on reaction cross 
sections. The collision trajectories have been studied using computer- 
plotter techniques. The interaction potential consists of a hard-sphere 
term plus the attractive polarization and ion-dipole terms. Multiple 
reflections corresponding to formation of long-lived complexes occur off 
the hard-core barrier at small ion-molecule separations. Multiple- 
reflection times are 10 to 500 times corresponding values for single re- 
flections. Hindered rotation of the polar targets (suggested by time 
history plots) is studied by computer-made movies. This technique also 
offers a convenient visual means for observing the relative motion of the 
ion-molecule pair. 

INTRODUCTION 

Capture cross sections have been calculated for ion-polar molecule 
collisions; the interaction potential consists of classical ion-dipole 
plus ion-induced dipole (polarizability) terms. These capture cross 
sections set a rough upper limit tp the reaction cross section. The ion- 
dipole energy transfer and formation of ion-molecule collision complexes 
have been studied using computer-plotter techniques. Mutual orbiting of 
the ion-molecule pair and hindered rotation of the dipole have also been 
studied via computer-made motion pictures. The latter is of particular 
interest for ion-molecule reactions since preferred dipole orientation 
favors certain reactions. This paper reports on results obtained for 
ion-molecule systems of particular interest in mass spectrometry.. 

COLLISION TARGETS 

The polar targets studied were chosen mainly because they have been 
studied in the mass spectrometer. Their molecular constants also cover 
an interesting range of collision parameters. Ion-molecule reactions 
involving CHgCN have been studied by Moran and Hamill.l The good 
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agreement between theory and mass spectrometer experimental values have 
been reported . 2 Large cross sections <v650 A 2 for CH 3 OH targets were 
reported in refs. 3 and 4 for low ion energies. Negative ion-molecule 
cross sections reported for H^O targets are roughly 10 3 cm 2 .-* 


The dipole moments, polarizabilities and moments of inertia as- 
signed to the three targets are given in Table I. CH 3 CN is a symmetric 
top whereas CH 3 OH and H 2 O were approximated as symmetric tops. The two 
moments of inertia assigned for each target were the largest moment and 
the moment about the symmetry axis. This approximation should not be 
critical for CH 3 OH; however, for H 2 O the neglected moment is approximately 
equal to the principal moment. ^ Previous comparison of capture cross 
sections for linear molecules and symmetric tops suggest that the cross 
section is independent of target geometry. ^ 


ION MOLECULE INTERACTION 


The ion-molecule interaction consists of the ion-permanent dipole and 
ion-induced dipole (polarizability) terms. The first term = -pe cos y 
where p is the permanent dipole moment, e is the electronic charge, y 
is the ion-dipole orientation angle and r is the ion-molecule separation. 
The Langevin polarizability term is -ae 2 /2r^ where a is the 

molecular polarizability which is assumed spherically symmetric. The 
separate cross section contributions are the Langevin cross section 


= n(2ae 2 /e) 1 ^ 2 = 


itb 2 


a) 


and the ion-dipole cross section. 


°b 


= itpe/e 


( 2 ) 


Assuming favorable orientation of the dipole, i.e., y = 0 , the cross 
section 


°M 


a L + a D = 


(3) 


is somewhat of a maximum. In equations (1) to (3), € is the relative 
translational energy. The ion and molecule are reflected off a hard- 
core potential at r = r . This reflection simulates ion-molecule re- 
pulsion due to the interacting electronic clouds. 
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COMPUTER APPROACH 

Collision quantities such as ion velocity and dipole rotational energy 
are plotted using values obtained by numerical integration. The time history 
plots are done in the coordinate system where the polar molecule is fixed at 
the origin.® The differential equations of motion are integrated using the 
variable step-length Runge-Kutta scheme of Ref. 7. Solutions of the equations 
include the coordinates for movie models of the colliding partners. The 
initial conditions and successful steps (for variables and time derivatives) 
are stored in a plotting array for both computer-made time history plots and 
movies. After the integration is completed the plotting subroutines are 
called. Each plotting array contains N + 1 values where N is the number of 
successful integration steps (typically 300 to 1000 ). 

MOVIE MODEL 

The coordinate solutions are used to draw projections of the ion and 
polar molecule models in the center-of-mass system. In fig. 1, a sample 
movie frame shows models of the ion and dipole with a clock in the upper 
right hand corner. ^ Each frame is traced on a cathode -ray-tube and photo- 
graphed using the IBM 360/67 computer with the GDC DD280 plotter. A hard 
sphere barrier is located at an ion-molecule separation r c equal to the 
sum of the ion and molecule radii. The plotted radius of each particle in 
the model is varied proportionately to the particle distance above (expanded) 
and below (contracted) the X-Y plane. The clock measures real time and moves 
at different rates during a collision since a variable step-size integration 
routine is used. 


INITIAL CONDITIONS 

The initial conditions on the dependent variables and their first time 
derivatives are obtained using a random number generator. Trajectories were 
studied for ion velocities of 5x10 , 10^, and 2x10^ cm sec"*-. Target rotators 
were chosen from a heat bath at rotational temperature T^ = 500 K. The 
estimated reflection distances for the collision complex studies were 2, 2.5, 
and 3.0 A for CO, HC1, and CH 3 CN, respectively whereas 2 A was used as a 
capture distance in all cross section calculations. It should be noted that 
these values are somewhat (25 percent) less than the Van der Waals radii of 
the neutral species. 0 


RESULTS 

In all numerical calculations the energy equation and behavior of 
trajectories with time reversal were used as checks on the accuracy of the 
integration routine. 
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Capture Cross Sections 

For each impact parameter 36 to 50 trajectories were calculated for 
the capture cross section; these trajectories were truncated once r £ 2 A. 
Unlike Langevin collisions there is no all-or-nothing answer to the question 
of capture in ion-permanent dipole collisions for fixed initial impact param- 
eter. For fixed values of energy and impact parameter there is a probability 
that the system will arrive at an ion-molecule separation corresponding to 
capture. The fraction of collisions resulting in such a prescribed "minimum" 
separation we call the capture ratio. Such a fraction can be correlated with 
the ion-molecule reaction probability. 

CH 3 CN + HC1. - Results of the capture ratio C R versus b 2 are shown 

in Fig. 2 for the C^CN"*" + CH 3 CN collision at two different relative veloc- 
ities. Results are shown for HC1 + + HC1 at two different rotational energies 
and two ion velocities in Figs. 3(a) and (b). The plot of the Langevin "all- 
or-nothing" capture ratio is a step function in both Figs. 2 and 3. 


A reasonable definition of the "capture collision" cross section is 


<y c =n C R (b) a(b 2 ) 


where b Q is the impact parameter at which C R becomes zero . The integral 
is simply the area under the C R versus b 2 curve. The calculated value 
b 2 is compared with b R (eq. ( 1 )), b^ ; (eq. ( 2 )), and the sum (eq. (3)). 
The cross section results are relatively insensitive to the choice of cutoff 
ion-molecule separation. (equal to 2A for HCl and CH 3 CN) . Figures 2 and 3 
indicate that the shape of the C R plots is sensitive to both rotational 
energy and ion velocity. The very large effect of the dipole is clear from 
comparison with the Langevin prediction for a molecule with the polarizability 
of CH 3 CN. For HCl the effectiveness of the dipole is far less than the 
"adiabatic" or maximum effect. 

In Fig. 4(a) the calculated cross section for CH^CN is shown as a 
function of energy and compared with various analytical estimates. The 
most obvious feature of this figure is how much larger the calculated 
cross section is than the Langevin cross section. The calculation also 
shows that the effect of the dipole is not as great at low energies as 
predicted on the basis of the adiabatic approximation. The capture cross 
section as assumed by the experimentalists is the integral of a Q 

over all ion energies. This cross section is compared with the calculated 
capture cross section Q c and observed reaction cross section Q R in 
Fig. 4(b). 
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CH 3 OH. - The plots of Cr versus are given for the CH 3 OH (sym- 

metric top approximation) + CH 30 H + collision in Figs. 5(a) and (b) for ion 
velocities of 5x10^ and 10 ^ cm sec”^-. The corresponding cross section a Q 
is plotted versus relative translational energy in Fig. 5(c). It is clear 
that this upper limit to the reaction cross section is somewhat less than 
the reported values even at thermal energy (€^_ = kT fc for T t = 500 K) . 

The 0-H bond does not lie along the symmetry axis of CH 3 OH (i.e., CH 3 OH is 
not quite a symmetric top). However, comparison of capture cross sections 
for linear molecules and symmetric tops suggests that additional degrees of 
rotational freedom have little effect on the capture probability .' 7 The 
agreement between theory and experiment (refs. 1, 3, and 4) is not as 
satisfactory as for the CH 3 CN system but this may be due to systematic 
experimental error . 


H” + H 2 O. - Stockdale et al.^ report cross section values of 5x10"^ cm^ 

(lower limit) at ion energies of € ]_ = 0.2 eV. The maximum cross section 
0^ computed for H" + H 2 O is only 2.5x10”^ cm^ (250 A^) at = 0.2 eV 

(corresponding to an H" velocity of 2x10^ cm sec“^). However, the numerically 
calculated cross section is only it (36 A^) = 113 (stee Fig. 6). The ratio 
of experimental to calculated cross section is thus roughly a factor of 45. 
This discrepancy does not appear explainable on the basis of the usual ex- 
perimental errors. More careful experiments are clearly required to obtain 
reliable absolute cross section data. 


Collision Lifetimes and Hindered Rotation 

Trajectories of the ion-molecule pair have been calculated as well as 
capture cross sections.8>9,ll Trajectory results indicate that there is a 
rapid interchange of energy between translational and rotational energy. 

This energy transfer suggests the possibility of long-lived collision 
complexes. Experimental evidence points to the existence of long-lived ion- 
molecule collision complexes. 12,13 Numerical lifetimes are only 10” 1-0 sec. 

To explain the long experimental lifetimes (as long as 10”? sec) requires 
consideration of additional mechanisms including participation of internal 
degrees of freedom. Such long-lived complexes must, therefore, be formed 
by inelastic collisions. 

For trajectory studies it is necessary to treat the short-range inter- 
actions in some fashion. We have used the following extreme approximation; 
a hard inner core, i.e., at r = r c the potential is stepwise infinite. The 
inward and outward trajectories are symmetric (specular reflection) for the 
simple Langevin potential. 

The situation is different for an ion-molecule collision involving a 
molecule with a dipole moment. The rotating dipole alters the potential for 
the outward trajectory, outer turning points can be introduced and multiple 
reflections occur. For parameters of interest in actual collisions as many 
as 2000 reflections have been observed numerically with collision times t 
as large as 10 ^ times the single reflection period.®’^ 
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Results of Trajectory Calculations 

Results of computer calculations of ion-molecule collisions are shown 
in Figs. 7 to 15. These time history plots have been made using the IBM 
360/67 computer with the CDC DD280 plotter. Figures 7 and 8 show the vari- 
ation of velocity and rotational energy as a function of separation. The 
collision shown in Fig. 7, for HC1 + NC> 2 + , involves a single reflection. 

The velocity increases as the separation decreases to the critical value 
where the potential energy is a minimum. The rotational energy changes 
drastically: there is hindered rotation at the small separations. 

The trajectory of Fig. 8 , for CO + Ar + , involves multiple reflections. 
Variation of the outer turning point is shown in this figure in the rotational 
energy plot. 

The simplest systems on which our computations are based are essentially 
triatomic (i.e., diatomic rotors and monatomic ions). The trajectories do not 
lie in planes as is shown in Fig. 9 obtained from calculations for CO + Ar . 
This figure shows a multiple reflection collision. 

Symmetric top dipolar molecules have one more degree of freedom than 
rigid rotors and their collisions are more complicated. Their trajectories 
involve more multiple reflections than rigid rotors. Figures 10 and 11 show 
trajectories for collisions of CH 3 CN with its parent ion (i.e., as regards 
ion mass). These figures indicate a more effective transfer of energy between 
rotational and translational degrees of freedom than for the linear molecule. 

Figures 12 to 15 show variations of the polar and azimuthal angles for 
trajectories involving multiple reflections for HC1 + CH 3 CN + CHgCN*, 

and CO + Ar + . Figure 15 shows how distant some of the outer turning points 
in multiple reflection trajectories can be. One turning point for CH 3 CN + 
CH 3 CN + occurs at 22 % separation. So-called "sticky" collisions involving 
alkyl halide molecules have been studied in the mass spectrometer; these 
collisions have large cross sections in polar systems. ^ These collisions 
should be similar to the CH 3 CN - CE^ON"*" collisions. 

It should be noted that multiple reflections are also to be expected 
for collisions involving polyatomic molecules without dipole moments since 
actual short range potentials probably have sufficient angular dependence 
in most cases to produce in varying degree the same effect as dipoles. 

Lifetime Considerations 

Theories of the lifetime of complexes have been considered-^’ 16 an( j can 
be applied to these cases. These theories assume that the energy of the 
complex is completely randomized over certain participating degrees of free- 
dom, and the rather complicated nature of our calculated multiply-reflected 
trajectories would suggest that such an assumption is reasonable for our 
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cases. Application of any theory requires more knowledge of the potential 
functions of the complexes than we have at present and the discussion pre- 
sented here can be only fragmentary. The simplest "theoretical" expression 
for the lifetime of a complex with enough energy to decompose is^ 

t = t 0 [e/(E - Eh)] 3 * T 0 [(Eh + e)/e] S (5) 

where t 0 is a time parameter made up of physical constants of the complex; 

E is the energy of the complex with its zero point energy taken as the zero 
of energy; Eh is the magnitude of the binding energy of the complex; s is 
a parameter determined by the number of degrees of freedom of the system. 

We have that E = Eh + £ where € is the relative energy of colliding 
pair. Thus t 0 is a constant of the same magnitude as the collision time 
without reflections; i.e,, t 0 lO - ^ s . Eh is in the range of 0.5 to 
several eV. For the various collisions we are considering here, s varies 
from 0.5 to 1.0. Since (Eh + £)/ £ varies from about 10 to 500, we expect 
that the collision time should be increased by a factor of 10 to several 
hundred by the reflection phenomenon. 

Computer-Made Movies of Ion-Dipole Collisions 

The computer-plotter studies have been extended to the making of motion 
pictures of ion-dipole single- and multiple-reflection collisions. The 
collision movies provide instantaneous visual correlation of the relative 
translational motion of the pair and the dipole rotational motion. One of 
the most interesting phenomena observed in the movies is the hindering of 
the rotational motion of the dipole. The three polar rotors CO, HC1, and 
CH 3 CN are seen to be hindered by the incident ion to varying degrees. U 

Figure 16 is a superposition of movie frames of a CH 3 CN capture collision 
involving only one reflection. It provides a correlated history of ion-dipole 
interaction. The results are not presented at equal time intervals since a 
variable step-size is used. The post-reflection hindering of the CH^CN rotor 
at 5 to 15 A is demonstrated in the movie. It has been suggested that prefer- 
ential orientation of the negative end of the dipole toward the positive ion 
will favor a specific chemical reaction. 

Figure 17 was obtained in the same manner as Fig. 16 except it shows 40 
frames superimposed for a multiple reflection collision. This collision 
sequence does not show hindered rotation but presents an interesting history 
of translational mutual orbiting. This mutual orbiting of the ion-molecule 
pair is nearly symmetric for this six reflection collision. This particular 
collision results in approximately 0 ° scattering in the center-of-mass system. 
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Vibrational Effects 

Numerical studies of ion-dipole collisions have been extended to 
calculating the effects of vibrational degrees of freedom on both multiple 
reflection probability and average collision lifetime t«. The dipole 
moment of the polar molecule is assumed Gaussian about the equilibrium 
bond separation of the oscillator. Both time history plots of collision 
variables and computer-made motion pictures are used to study the inter- 
action of the ion with the oscillating dipole. Initial results have been 
obtained chiefly for loose (thermal) oscillators with energies of 0.036 eV. 

The results indicate such oscillators actually lower the fraction of multiple 
reflection collisions for CO but increase it for CH 3 CN at impact parameters 
greater than 8 5L The collision lifetimes are somewhat insensitive to the 
presence of the oscillator for the 300 cases studied. Stiff oscillators 
(50 cases studied) generally lower the value by cooling; this result 

is unphysical for targets vibrating initially in the ground state. 

CONCLUDING REMARKS 

Numerical calculations have been made on several different features of 
ion-polar-molecule collisions. Calculated cross sections are in good to 
satisfactory agreement with experiment for CHgCN and CHgOH targets. However, 
there is considerable disparity between predicted and observed values for 
negative ion-^O reactions. There is a clear need for more reliable absolute 
cross section data; the numerical results also suggest that the rotational 
temperature dependence of the cross section should be investigated. 

The existence of long-lived ion-molecule complexes ( t » ICT^ - 10 “-^ sec) 
has been demonstrated for a rotating polar target with and without vibrational 
effects. A more sophisticated model incorporating additional internal degrees 
should predict collision complex lifetimes in better agreement with experiment 
('vlCT 7 sec). (This conclusion is based on simple lifetime predictions.) 
Preliminary calculations which take thermal harmonic oscillators into account 
remain somewhat inconclusive with regard to lifetime effects. 

Time history plots of ion-dipole orientation angle have suggested that 
considerable hindering of the dipole occurs. Gomputer-made movies have proven 
valuable in demonstrating hindered rotation in ion-dipole collisions via 
instantaneous correlation of translational and rotational motion. Super- 
imposed movie frames provide time histories of this hindered rotation and 
mutual orbiting of collision partners. An ultimate objective of this work is 
the application to more complicated processes such as ion-molecule reactions 
and charge exchange. In such processes it is quite clear that relatively 
long-lived complexes must be involved. 
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TABLE I 


Polar 

target, 

molecule 

Dipole a 

moment, 

D.U. 

Po la ri zabili tv - 
p 3 
A 

Principal 
moment of inertia, 
l2_( gm-cm 2 )xl0^° 

Ratio of 
moments of 
inertia, 

12A1 

CH3CN 

3.92 

3.8 

9.12X10 1 

0.06 

CH3OH 

1.67 

2.4 

3.58X10 1 

.15 

h 2 o , 

1.86 

2.5 

1.93 

1.51 


a l Debye unit (D.U. ) = 10 “^® esu -cm. 
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Figure 2. - Variation of capture ratio with square impact parameter for 
CH 3 CN + + CH 3 CN system. The symmetric top target rotators are dis- 
tributed at temperature T r = 500 K. Bars indicate variation of calcula- 
tion with different random number sets. 



CAPTURE RATIO, C. 


HC! + + HC1 



(a) v = 10 5 cm/sec. <b) v □ 5. 5X10 4 cm/sec. 


Figure 3. - Variation of capture ratio with square of impact parameter for linear HCI target 
rotators distributed at temperatures T R * 500 and 1000 K; 50 collisions were studied per 
point. 



GROSS SECTION, INTEGRATED OVER MONOENERGETIC CROSS SECTION, o, A 2 



(a) Comparison of numerically calculated cross sections a c for methyl cyanide- 
parent ion collision with various theoretical cross sections. All cross sections 
are plotted as functions of ion translational energy for one rotational tempera- 
ture T r = 500 K ; electronic polarizability a, 3. 8 A 3 ; dipole moment p, 3.92 
Debye units. 



MAXIMUM ION ENERGY, e m , eV 


(b) Comparison of numerical capture cross section Q c with 
experimentally assumed capture cross section Q max and 
observed reaction cross section for methyl cyanide-parent 
ion collision as function of maximum ion energy € m . 


Figure 4 




CAPTURE RATIO, C 




Figure 5. - Variation of capture ratio with free-particle impact parameter for rotational temperature 
of 500 K. Methyl alcohol parent-ion collision. 



ION-MOLECULE RELATIVE ENERGY, <■, eV 

(c) Comparison of integrated numerical 
capture cross section Q c with exper- 
imentally observed capture cross 
section Q R for methyl alcohol 
parent-ion collision. 

Figure 5. - Concluded. 



ROTATIONAL ENERGY VELOCITY SCALED CAPTURE RATIO C 

SCALED TO kT R = 0.026 eV TO INITIAL VALUE OF 

5.5X10 4 CM SEC" 1 k, In 



Figure 6. - Variation of capture ratio with free-particle 
impact parameter for rotational temperature of 500 K. 
FT + H 2 0 capture collision; « * 0. 2 eV ; 
v * 6.36xl0 5 cm sec" 1 ; 50transectories per point. 



ION-MOLECULE SEPARATION IN A 


Figure 7. - Variation of ion velocity and polar molecule 
rotational energy during N0 2 + HC! single reflection 
capture collision. 



Figure 8. - Variation of ion velocity and polar molecule 
rotational energy during Ar + + CO multiple reflection 
capture collision. 


Figure 9. - Variations of polar angle 6 for translational 
motion of Ar + relative to CO molecule during multiple 
reflection capture collision. 









Figure 10. - Variation of ion velocity and polar 
molecule rotational energy during methyl 
cyanide-parent ion multiple reflection 
capture collision. 



Z AXIS ION-MOLECULE 

SEPARATION IN A 

Figure 11. - Variations of dipole moment vector and ion- 
dipole orientation angle during C^CN-CF^CN* multiple 
reflection capture collision with several turning points. 



Figure 12. - Variation of polar angle for translational motion 
of CH 3 CN relative to CH 3 CN molecule during multiple re- 
flection capture collision. 
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ION-MOLECULE SEPARATION IN A 

Figure 13. - Variations of aximuthal angle (p for transla- 
tion motion of NOj relative to HCI molecule during 
multiple reflection capture collision. 












Z AXIS 

Figure 14. - Variation of ion projections tracing transla- 
tional motion in Ar + + CO capture collision with 
multiple reflections. 



Figure 15. - Variation of azimuthal angle (p for translational 
motion of CH 3 CN* relative to CH 3 CN molecule during 
multiple reflection capture collision with maximum turning 
point of 16 A. 



Figure 16. - Superposition of 15 movie frames 
for single reflection CH 3 CN parent-ion collision, 
ion-molecule pair are designated in pre-reflection 
positions; motion is traced during and after re- 
flection. 
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Figure 17. - Superposition of 40 movie frames 
demonstrating mutual spiraling of ion and 
molecule in CH^CN-parent ion multiple 
reflection collision with zero degrees net 
scattering. 
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